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Abstract. We investigated the relation between the centroid frequency of the quasi-periodic oscillation observed 
' in the power density spectra of a sample of galactic black-hole candidates with the power-law photon index 

t-H , obtained from spectral fits. Our aim is to avoid inner accretion disk radius determination directly from spectral 

fits, given the uncertainties of the absolute values obtained in that way, but to base our analysis on the likely 
association of QPO frequency to a characteristic radius. We used archival RXTE data of GRS 1915+105 and 
published parameters for GRO 1655-40, XTE J1550-564, XTE J1748-288 and 4U 1630-47. While for low values 
of the QPO frequency, the two parameters are clearly correlated for each source, there is evidence for a turnoff in 
the correlation above a characteristic frequency, different for different sources. We discuss the possible nature of 
f*^ . this turnoff. 
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■ 1. Introduction also Borozdin et al. 1999). The numerous observations of 

' GRS 1915+105 made with the instruments onboard the 

k1 " In 1992 ' th<3 WATCH all " sky monitor onboard GRANAT Rossi x _ my Timing Explorer (RXTE) satellite show an 

• discovered the X-ray source GRS 1915+105 as a extremely complex variability in the X-ray band, which 

g ! bright transient (Castro-Tirado et al. 1992). Mirabel & Belloni at aL (1997a, 1997b) interpreted as an oscilla- 

Rodrfguez (1994) reported the observation of radio rel- tion of thc inncr rcgion of the accret ion disk caused by 

ativistic expansions interpreted in terms of collimated a thermal-viscous instability. Belloni et al. (2000) intro- 

emissions of matter, which made GRS 1915+105 the first duced a model-independent classification of the different 

galactic object with superluminal jets, for which the term types of x _ ray variab iii ty m three basic spectral states, 

"microquasar" was coined. This source was originally sup- Similarly to more conventional Black Hole Candidates 

posed to host a black hole because of its high X-ray lu- (BHC), the energy spectrum of GRS 1915+105 can be fit- 

minosity, and for its similarity with another superluminal te d with the superposition of a soft thermal component, 

source, GRO J1655-40, for which it was possible to ob- mode i ed as a d i s k-blackbody (Mitsuda et al. 1984), plus 

tain a dynamical estimate of the mass (M=6.3±0.5M Q ) a harder C0mp0Ilcn t, modeled as a power law with a high 

of the compact object (Bailyn et al. 1995; Greene et al. energy cutoff ( Ta naka & Lewin 1995). The state transition 

2001). No optical counterpart of GRS 1915+105 was found can be attributed to changes in the relative contributions 

because of the high Galactic extinction and only an IR of these variable components, again similarly to what ob- 

counterpart was detected (Mirabel et al. 1994). Recently, served in other BHCs (Belloni 1998) 
an orbital period of 33.5 days was discovered with IR spec- 
troscopic observations, leading to a dynamical estimate of Detailed timing analysis of the RXTE data (Morgan 

14±4 Mq for the compact object (Greiner et al. 2001, see et al _ 1997 ) revealed three types of Quasi Periodic 

Oscillations (QPO): a 1-10 Hz QPO with variable cen- 

Send offprint requests to: T. Belloni troid frequency similar to that observed in other BHCs 
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(see van der Klis 1995), two quasi-stable QPOs with cen- 
troid frequencies ~40 and ^67 Hz (Morgan et al. 1997; 
Belloni et al. 2001; Strohmayer et al. 2001a), and variable 
low frequency QPOs (10~ 3 — 1 Hz), which can be identified 
with transitions between the states of Belloni et al. (2000). 
The presence of the 1-10 Hz is always associated with a 
spectral state C (Belloni et al. 2000). The 1-10 Hz QPO 
centroid frequency is well correlated with both disk and 
power-law spectral components (Muno et al. 1999, Reig et 
al. 2000, Rodriguez et al. 2002) and with the X-ray flux 
(Chen et al. 1997; Markwardt et al. 1999). 

Several theoretical models associate the origin of this 
QPO with the innermost region of the disk (Miller et 
al. 1998; Nowak & Wagoner 1991; Titarchuk et al. 1998; 
Stella & Vietri 1998; Stella et al. 1999; Varniere et al. 
2002, Psaltis & Norman 2002). Di Matteo & Psaltis (1999) 
made use of the fact that the fastest variability timescale 
at any radius around a compact object is the Keplerian 
orbital frequency to obtain upper limits on the radius at 
which the 1-10 Hz QPO is produced. As it is apparent 
in their work, there is a general positive correlation, con- 
sidering different systems, between the photon index of 
the power-law component in the energy spectrum and the 
centroid frequency of the QPO. Sobczak et al. (2000a) 
in their analysis of RXTE data of two transient BHCs, 
XTE J1550-564 and GRO J1655-40, examined the cor- 
relation between the spectral parameters, inner disk ra- 
dius and power-law photon index, with the QPO centroid 
frequency. Surprisingly, they find correlations of opposite 
sign for the two sources. Rodriguez et al. (2002) analyzed 
the correlation between QPO centroid frequency and inner 
disk radius (determined from X-ray spectral fits). 

In this paper, we concentrate on the correlation be- 
tween power-law photon index and QPO centroid fre- 
quency. We present the results of timing and spec- 
tral analysis of a number of RXTE observations of 
GRS 1915+105, focusing on this correlation. We also 
consider published data relative to four other systems: 
GRO J1655-40 (Remillard et al. 1999a; Sobczak et al. 
2000a), XTE J1550-564 (Sobczak et al. 2000a,b; Homan 
et al. 2001), XTE J1748-288 (Revnivtsev et al. 2000), 
4U 1630-47 (Trudolyubov at al. 2001; Tomsick & Kaaret 
2000). We then compare the results in term of possible 
theoretical models. 



2. Data analysis 

Our study of GRS 1915+105 is divided in two parts, ac- 
cording to the different classes of observations that we con- 
sidered. In the first part, we analyze the so-called plateau 
observations (Fender et al. 1999; Reig at al. 2000), which 
consist in long intervals of persistent state C (class x m 
Belloni et al. 2000). The second part is dedicated to obser- 
vations that show large variability in the X-ray flux and 
at least one state-C interval. Again following the classi- 
fication of Belloni et al. (2000), we selected observations 
corresponding to classes a, /3, and v. 



2.1. Plateau observations 

Plateau intervals show persistent QPOs that are very eas- 
ily detectable in the Power Density Spectrum (PDS). We 
based our choice of appropriate observations upon the ex- 
tensive analysis of Reig et al. (2000), which allowed us to 
identify the larger possible range of centroid frequencies of 
the QPO. We selected 7 RXTE/PCA observations of GRS 
1915+105 (for a total of 15 satellite orbits), corresponding 
to class x m the classification of Belloni et al. (2000). The 
observation log is reported in Table [|. 

The plateau observations are relatively straightforward 
to analyze from the point of view of timing analysis, be- 
cause of the stable positive correlation between the QPO 
frequency and count rate (see Reig et al. 2000). Therefore, 
no special timing selection is necessary and it is possible 
to accumulate all the data to produce one PDS and one 
energy spectrum. To produce the PDS, for each observa- 
tion we used high-time resolution data in the 2-13 keV 
range. Some of the observations have a time resolution of 
1/512 s, limited to 1/128 s for the others (see Table @). 
We divided each observation interval (orbit of the satel- 
lite) in intervals 16 s long and produced an average PDS 
per observation summing the PDS of the individual seg- 
ments. We subtracted the contribution due to Poissonian 
statistics and the Very Large Event window (Zhang et 
al. 1995, Zhang 1995) from each PDS and renormalizcd 
them to squared fractional rms per Hertz (see Belloni & 
Hasinger 1990). Every PDS was then logarithmically re- 
binned and fitted to find the centroid frequency of the 1-10 
Hz QPO, which is very sharp and often shows higher har- 
monics and even a sub-harmonic component (as shown in 
Fig. [ljfor three different observations). A model consisting 
of the superposition of a number of Lorentzians (from 5 
to 9 depending on the observation) was used to fit every 
feature present in the power spectrum (see Belloni et al. 
2002). Because of the relative strength of the QPO peak, 
the centroid frequency did not change appreciably when 
using different models for the continuum noise: the QPO 
frequencies we found with this analysis are in close agree- 
ment with those from Reig at al. (2000). The resulting 
QPO centroid frequencies are listed in Table [l]. 

In Class x observations, the source is in a long steady 
state C (Belloni et al. 2000). As shown by Belloni et 
al. (2000), this corresponds to an undetectable thermal 
component in the PCA energy spectrum. For each ob- 
servation interval we extracted PCA energy spectra in 
the range 3-25 keV, averaged over the whole interval, 
from spectral (Standard2) data. For each spectrum, we 
subtracted the background estimated with pcabackest 
v. 2.1e, and created an appropriate detector response 
matrix with pcarsp v. 2.43. A systematic error of 1% 
was added to the spectra to account for remaining un- 
certainties in the detector calibration. Our starting model 
for the spectral fits consisted of a power law, corrected for 
interstellar absorption, plus a Gaussian emission line to 
take into account an excess at 6.4 keV. With this model, 
it was not possible to lower the reduced x 2 below 2. We 
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Obs.N Observation ID MJD PCA(c/s) AT v QPO Tim.Res. Y N H E c xl 



1 


10408 


- 01 


- 23 


- 00a 


50278.4921 


9735 


3072 


3.491 ± 0.005 


1/512 


2.62 ± 0.01 


6.7 


24.7 


1.7 


2 


10408 


- 01 


- 23 


- 006 


50278.5585 


9710 


3540 


3.617 ± 0.004 


1/512 


2.66 ± 0.01 


6.9 


25.4 


1.7 


3 


10408 


- 01 


- 23 


- 00c 


50278.6266 


10458 


3401 


4.199 ± 0.009 


1/512 


2.76 ± 0.01 


7.1 


28.0 


1.7 


4 


10408 


- 01 


- 24 


- 00a 


50280.1702 


8924 


2497 


2.228 ± 0.008 


1/512 


2.38 ± 0.04 


6.4 


25.1 


1.7 


5 


10408 


- 01 


- 24 


- 006 


50280.2266 


8901 


3479 


2.293 ± 0.007 


1/512 


2.41 ± 0.04 


6.3 


23.0 


1.6 


6 


10408 


- 01 


- 24 


- 00c 


50280.2933 


8956 


2998 


2.534 ± 0.008 


1/512 


2.46 ± 0.04 


6.6 


23.0 


1.4 


7 


10408 


-01 


- 27 


- 00a 


50290.5774 


7887 


2421 


0.642 ± 0.004 


1/512 


1.88 ± 0.04 


6.9 


16.2 


1.3 


8 


10408 


-01 


- 27 


- 006 


50290.6322 


7926 


3480 


0.628 ± 0.004 


1/512 


1.82 ± 0.01 


6.7 


15.1 


1.2 


9 


20402 


-01 


- 49 


- 01 


50730.3949 


8198 


3463 


2.637 ± 0.005 


1/512 


2.45 ± 0.01 


5.4 


29.7 


1.8 


10 


20402 


-01 


- 50 


-00 


50735.5474 


6579 


1962 


0.835 ±0.010 


1/128 


1.92 ±0.02 


5.7 


21.5 


1.9 


11 


20402 


-01 


- 50 


- Ola 


50737.4047 


6608 


3007 


1.014 ±0.006 


1/128 


1.98 ±0.01 


5.7 


22.1 


1.6 


12 


20402 


-01 


- 50 


- 016 


50737.4766 


6373 


2620 


1.078 ± 0.005 


1/128 


2.00 ±0.01 


5.5 


20.0 


1.6 


13 


20402 


-01 


- 52 


- 00a 


50746.5509 


6657 


2207 


1.405 ±0.007 


1/128 


2.09 ± 0.04 


5.5 


25.0 


1.8 


14 


20402 


-01 


- 52 


-006 


50746.6190 


6668 


2079 


1.468 ± 0.006 


1/128 


2.10 ± 0.04 


5.5 


25.1 


1.4 


15 


20402 


-01 


- 52 


-00c 


50746.6898 


6768 


1725 


1.609 ±0.007 


1/128 


2.16 ± 0.04 


5.6 


26.9 


1.7 



Table 1. List of plateau observations analyzed and results for timing and spectral parameters. Letters a, b, c correspond 
to different orbits in the same observation. The fifth column refers to the state-C length (in seconds) that we selected 
for our study. QPO frequencies are in Hertz, Nh X 10 22 cm -2 . Typical errors are: for Nh±0.1 and for E C ±0A5 



modified the model, adding a high-energy cutoff to the 
power law. The resulting spectral parameters are shown 
in Table 1. This model is consistent with that used by 
Trudolyubov et al. (1999) and Muno et al. (1999). Notice 
that Trudolyubov et al. (1999) found rather high cut- 
off energies (-100 keV), while Muno et al. (1999) found 
also much lower values, even below 10 keV. Our values 
are in the range 20-30 keV. Only 2 among the selected 
plateau observations (10408-01-27-00 and 20402-01-52-00) 
were also analyzed by Muno et al. (1999). They used an 
often-used model for Black Hole Candidates, which con- 
sists of a superposition of a multicolor disk-blackbody 
(thermal component) and a power law (without high en- 
ergy cutoff). Their spectral analysis results in a different 
power law photon index, ~ 0.7 grater than the values we 
found. However, they derive an extremely high inner disk 
temperature (~ 4 keV) and an inner disk radius consistent 
with zero (~ 2±2 km). These parameters do not represent 
a realistic physical picture of the source and therefore sug- 
gest that the inclusion of a disk-blackbody component is 
not physically necessary. Our choice of a simple power law 
with a high-energy cutoff is motivated by two reasons: the 
need of a simple model for efficiently performing a large 
number of fits, and the need to use a common model that 
allows comparison with values in the literature. In order 
to check that for high values of the spectral index T this 
model yields meaningful values, we fitted our spectra also 
with the bmc model in XSPEC, which has been used by 
other authors to describe the energy spectra of Black Hole 
Candidates in the framework of Comptonization models. 
These fits give values of T consistent with those from a 
simple power-law model, with a slight systematic excess 
of 3%. From Table |, one can see that there are variations 
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0.5 1 2 5 10 

Fig. 2. Plot of Power Law Photon Index versus QPO 
centroid frequency for the plateau observations of 
GRS 1915±105 studied. A positive correlation spanning 
nearly an order of magnitude is evident. Representative 
error bars are plotted for three points. 

in the best fit values of the column density (Nh) (see also 
Belloni et al. 2000; Klein- Wolt et al. 2002). 

In Fig. ^|we plot the power-law photon index versus the 
QPO centroid frequency. A positive correlation is evident. 

2.2. Class a, f3 and v observations 

In order to examine the correlation between QPO centroid 
frequency and power-law photon index, and in particular 
to extend it to higher values of the QPO frequency, we 
also analyzed observations in which GRS 1915±105 shows 
large variability on time scales longer than 1 second. While 
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belonging to different classes, they all have in common the 
presence of at least one state-C interval of more than 700 
seconds (see Belloni et al. 2000). The log of selected ob- 
servations is reported in Table ^. After a close inspection 
of the light curves produced from Standardl data (cor- 
responding to the full PCA energy band 2-60 keV), we 
identified the start and end times of the state-C inter- 
vals to analyze. Flux values during these periods typically 
range from 3000 to 40000 cts/s. For each observation, we 
checked the choice of intervals by producing a dynamical 
PDS (a series of PDS as a function of time), with which we 
verified that a 1-10 Hz QPO was present throughout the 
whole interval. The dynamical PDS also confirmed that 
the X-ray intensity of the source is strongly correlated 
with the QPO frequency (Chen et al. 1997, Markwardt et 
al. 1999, Reig et al. 2000). Therefore, following the PCA 
count rate, the QPO frequency first decreases rather fast, 
then starts increasing at a lower rate (see Fig. [|). 

To perform timing analysis, we inspected the same 
type of data (high-time resolution and low energy range) 
as for plateau observations. In this case however, we di- 
vided every state-C interval in segments with a length of 
4 s rather than 16, in order to better follow the rapid 
variations of the QPO centroid frequency, and produced a 



Obs.N Obs.ID MJD Class AT 



15 


10408 


- 01 


-44- 


00a 


50381.4942 


7/ 


1248 


16 


10408 


-01 


-44- 


00b 


50381.5655 


V 


1488 


17 


20186 


- 03 


-03- 


016 


50674.2347 




704 


18 


20187 


- 02 


-01- 


00b 


50575.9322 


a 


928 


19 


40702 


- 01 


-02- 


00b 


51320.9322 


V 


1120 



Table 2. List of the observations from classes a, (5 and 
v analyzed in this work. Letters a, b, c correspond to dif- 
ferent observation intervals within the same observation. 
The last column refers to the length (AT) of the state-C 
interval that we selected for our study 



PDS for every segment. After subtracting the contribution 
from the Poissonian statistics and the Very Large Event 
window, we summed these PDS in sets of four, rebinncd 
them logarithmically and fitted the result to find the QPO 
centroid frequency. As in the previous case, a model con- 
sisting of the superposition of a number of Lorentzians was 
used (see Belloni et al. 2002). We performed the summing 
process because we wanted to correlate timing and spec- 
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Fig. 3. Examples of power density spectra (left panels) and energy spectra (right panels) for 40702-01-02-00 (class v) 
corresponding to different time intervals. Start and end times for each panel are from the observation start time (Table 
2). 



tral results (which have a minimum time bin of 16 s), and 
this is possible only if we analyze intervals with temporal 
coincidence. 

Spectral parameters were obtained extracting energy 
spectra from Standard2 data, with an integration time of 
16 seconds, subtracting the background and creating the 
detector response matrix as for the plateau observations. 
We fitted every spectrum using the "standard" model for 
BHCs: the superposition of a multicolor disk-blackbody 
and a power law for the high energy tail. We corrected 
for the interstellar absorption, fixing Nh = 7 x 10 22 cm~ 2 
(Klein- Wolt et al. 2002) and added a Gaussian emission 
line with central energy fixed to 6.4 keV. With respect to 
the model adopted for the plateau observation, we needed 
to include a disk-blackbody (see Belloni et al. 1997b, 2000) 
and we had to fix Nh due to the limited statistics of the 
16s spectra. Again, we added in quadrature a systematic 
error of 1%. The fits give acceptable results, with a re- 
duced x 2 usually around 1. Therefore, we can trace the 
evolution (with a time step of 16 s) of the spectral param- 
eters during the state-C intervals. We can then compare 
the spectral and timing evolution by correlating the re- 
sulting parameters. As in the case of plateau observations 



we fitted selected spectra with the bmc model (see Sect. 
2.2), finding consistent results for T. 

Examples of both PDSs and energy spectra (with their 
best fit models) can be seen in Fig. y. The correlation be- 
tween QPO centroid frequency and power-law index for 
classes (3 and v is plotted in Fig. ^ (upper panel). The be- 
haviour is similar to \ observations, although there is a 
spread due to the lower statistics available from spectra 
averaged over only 16 seconds. At high QPO frequencies 
it is possible to notice a "turnoff" in the correlation for 
observations 15 and 16, while observation 19 (diamonds) 
does not show this feature. Figure (lower panel) shows 
the same correlation from classes a and v. While v data 
are compatible with those in the upper panel, a data show 
a correlation similar, but offset to lower photon-index val- 
ues. 



3. Other sources 

In order to compare our results for GRS 1915+105 with 
other systems, we searched the literature for similar data 
from other black hole candidates. We identified four sys- 
tems for which values for the QPO centroid frequency and 
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Fig. 4. Light curve and dynamical PDS for the z/-class 
observation 40702-01-02-00 (second interval). It is evident 
that the 1-10 Hz QPO, when present, is strictly positively 
correlated with counts (see text). In the PDS frame, fre- 
quency is on the y-axis while the power scale is shown in 
the chromatic column on the right. 

the power-law index are available, and for which there 
are at least a few different values. Notice that, as men- 
tioned above, we do not consider parameters like inner 
disk radius, which are difficult to compare between dif- 
ferent sources. Providing the adopted spectral models are 
compatible, the values of T are directly comparable, and 
so of course are the QPO centroid frequencies. 

3.1. GRO J 1655-40 

This source is a BHC with a dynamical estimate for the 
mass of the black hole of ~ 6.3M© (Bailyn et al. 1995; 
Greene et al. 2001). The system, a bright X-ray transient 
which appeared for the first time in 1994 and showed 
a subsequent outburst in 1996/1997, showed superlumi- 
nal radio jets (Hjellming & Rupen 1995). GRO J1655- 
40 is therefore another of the few objects belonging to 
the microquasar class. This object shows quite a few dif- 
ferent QPO peaks in a range between 0.1 and 450 Hz, 
(see Remillard et al. 1999a; Strohmaycr 2001, Rcmillard 
et al. 2002b). Among them, we are interested in the vari- 
able 14-22 Hz QPO (see Sobczak et al. 2000a). We con- 
sider here PCA/RXTE data from an intensive observa- 
tional campaign of the source during its second outburst 
in 1996-1997. The timing parameters (i.e., the QPO cen- 
troid frequencies) are taken from Table 3 of Remillard 
et al. (1999a), and the spectral parameters (i.e., the 
power law photon indexes) from Table I of Sobczak et 
al. (1999a). The energy spectra in the 2.5-20 keV energy 
range were fitted using a model consisting of a multicolor 




0.5 1 3 5 10 

QPO (Hz) 
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QPO (Hz) 

Fig. 5. Upper panel: power-law photon index vs. QPO cen- 
troid frequency for the observations of class and v stud- 
ied. Values for the plateau observations (from Fig. 2) are 
plotted for comparison (stars). Different symbols are used 
to identify different observations. [Triangles= obs. 15,16 
- Diamonds^ obs. 19 - Crosses^ obs. 17]. Lower panel: 
power-law photon index vs. QPO centroid frequency for 
the observations of class a and v studied. Different sym- 
bols identify different observations. [Full circles^ obs. 
15,16 - Empty circles^ obs. 19 - Squares^ obs. 18]. In 
both panels error bars at 90% confidence for three data 
points are shown as example. 

disk-blackbody plus a power law; a detailed description 
of the spectral analysis can be found in Sobczak et al. 
(1999a). The power density spectrum for each observa- 
tion was analyzed in the 2-30 keV energy range, obtain- 
ing the QPO central frequency fitting with a Lorentzian 
function plus a power law for the local continuum (see 
Rcmillard et al. 1999a). The plot of power-law index ver- 
sus QPO centroid frequency for GRO J1655-40 is shown in 
Fig. ^ (squares). There is a clear anti-correlation between 
these two parameters, as already shown by Sobczak et al. 
(1999a). 

3.2. XTE J1550-564 

This BHC is the brightest X-ray transient observed with 
RXTE (the 1998-1999 flare reached 6.8 Crab in the 2- 
10 keV band). In the PDS, it shows both low-frequency 
(0.8-18 Hz) and high-frequency (180-280 Hz) QPOs (Cui 
et al. 1999; Remillard et al. 1999b; Homan et al. 2001; 
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Remillard et al. 2002a, b). Wc consider here the low- 
frequency one (see Remillard et al. 2002a). The ob- 
servations presented here are the whole set made with 
PCA/RXTE before 1999 March 22 (Gain Epoch 3). The 
power-law photon indexes obtained from the spectral anal- 
ysis of these observations are taken from Table 1 of 
Sobczak et al. (1999b), while the QPO centroid frequen- 
cies from the PDS analysis are listed in Table 1 of Sobczak 
et al. (2000a). The spectral model used by Sobczak et al. 
(1999b) to fit the 2.5-20 keV energy spectra consists again 
of the superposition of a multicolor disk-blackbody and 
a power law, with corrections for interstellar absorption 
and a Gaussian emission line at ~ 6 — 7 keV. The centroid 
frequency of the QPO was found fitting the peak with a 
Lorentzian function and the continuum with a power law. 



& * 



the smooth decrease in flux until the first small flare, a 
'flat' interval where the flux remains rather stable, with 
small flares, and a 'decay' interval where the flux decays 
at the end of the first part of the outburst. After the last 
observation examined here, no QPO was detected for a 
long time, until the flux rose again to the second part of 
the outburst (see Homan et al. 2001). 

In Fig. 8 the behaviour of XTE J1550-564 in the QPO- 
r plane is shown with different symbols for the four differ- 
ent intervals. One can see that the first two intervals trace 
a very smooth curve in Fig. 8, the second interval trac- 
ing back the path of the first. The scattered points corre- 
spond to the remaining two intervals, and it is interesting 
to note that all these points are located below the others. 
What is different between the first two and the second two 
classes of observations that have a different behaviour in 
this plot? From Sobczak et al. (2000b) one can see that 
a clear dividing line between the two groups is the ratio 
of the power-law flux over the total flux, estimated from 
the observed PCA spectra. The first observations have a 
lower value of the ratio. Also notice that the points that 
contribute to the scatter in Fig. 8 are the ones from the 
'decay' interval, which contains QPOs classified as A/B in 
Homan et al. 2001, while the QPO observed in the pre- 
vious intervals is classified as C (notice that these QPO 
classes have nothing to do with the A/B/C classes defined 
in Belloni et al. 2000). 



Fig. 6. Power-law photon index vs. QPO centroid fre- 
quency for XTE J1550-564 (circles) and GRO J1655-40 
(squares) . While for the first source there is a general pos- 
itive correlation between the two parameters, for the sec- 
ond one the correlation is negative. However, one can see 
that data points of GRO J 1655- 40 seem to extend the 
"turnoff" curves of XTE 31550-564 (see text). 

Data obtained from the literature for XTE J 1550- 
564 are plotted in Fig. || (circles). They show a pos- 
itive correlation between power law photon index and 
QPO centroid frequency very similar to the behaviour of 
GRS 1915+105 (see Fig. | and |). However, while the 
correlation is very good at low frequency, at higher QPO 
values there is a larger spread. Surprisingly, the points 
from GRO J1655-40, although showing a negative corre- 
lation, are placed on the natural continuation of the upper 
points for XTE J1550-564. As one can see from Cui et al. 
(1999), Remillard et al. (2002b) and Homan et al. (2001), 
the low-frequency QPO in XTE J1550-564 does not main- 
tain the same properties throughout the whole outburst. 
It is therefore natural to examine its behaviour in separate 
intervals of the light curve. 

We subdivide the light curve in four separate sections 
(see Fig. 7): a 'rise' interval from the first observation up 
to the peak flux, a 'down' interval from the peak following 



3.3. XTE J 1748-288 

This BHC shows superluminal radio jets and belongs to 
the class of galactic microquasars, for the same reasons 
as GRS 1915+105 and GRO J1655-40 (Smith et al. 1998; 
Hjellming et al. 1998; Revnivtsev et al. 2000; Kotani et 
al. 2000; Miller et al. 2001). During the 1998 outburst the 
source was observed to display spectral and timing prop- 
erties typical for BHCs systems (Revnivtsev, et al. 2000). 
A QPO feature was found in the PDS at 0.5 Hz and 20- 
30 Hz. Here we consider spectral and timing parameters 
given in Table 2 and 3 of Revnivtsev et al. (2000) who re- 
ported an analysis of the PCA/RXTE observations during 
the 1998 outburst. They fitted the energy spectra between 
3 and 25 keV using the standard model for BHCs con- 
sisting of a multicolor disk-blackbody plus a power law, 
with low energy interstellar absorption and a Gaussian 
emission line at ~ 6.5 keV. For the seven observations 
considered here, Revnivtsev et al. (2000) fitted the PDS 
using the sum of a flat-topped band limited noise compo- 
nent, a power law component, and a Lorentzian function 
to describe the QPOs. The resulting points in the QPO-r 
plot arc shown as empty diamonds in Fig. 9. There is a 
clear positive correlation between T and QPO, but it is 
not identical to the same correlations of the three sources 
previously discussed. Notice that the QPO frequencies are 
the highest considered here. 
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Fig. 7. Bottom panel: light curve for XTE J1550-564 
as derived from Table 1 of Sobczak et al. (2000b). 
Observations span from 7th September 1998 to 20th 
November 1998. Different symbols correspond to the 
different intervals marked in the panel (see also Fig. 
8). Middle panel: corresponding QPO frequencies (from 
Sobczak et al. 2000a). Top panel: corresponding power- 
law slopes (from Sobczak et al. 1999b). 



3.4. 4U 1630-47 

This X-ray transient is considered a BHC for which the 
optical thin radio emission, observed during the 1998 out- 
burst, suggests the presence of radio lobes (Hjellming et 
al. 1999). We obtained timing (QPO centroid frequency) 
and spectral (power-law index) parameters of the obser- 
vations made during the decay of the 1998 outburst from 
Table 4 and 5 of Tomsick & Kaaret (2000). The timing 
analysis was made in the 2-21 keV energy band fitting 
the PDS with a fiat-top component for the continuum and 
a Lorentzian function for the QPO feature. The centroid 
frequency was between 0.2 and 3.4 Hz. The energy spec- 
tra were fitted in the 2.5-20 keV range using a power law, 
a multicolor disk-blackbody soft component, a Gaussian 
emission line at ~ 7 keV, the interstellar absorption and a 
broad iron absorption edge. For the peak of the outburst, 
we used the results of Trudolyubov et al. (2001). Here tim- 
ing analysis was performed from data in the 2-13 keV en- 



Fig. 8. Plot of Power Law Photon Index versus QPO cen- 
troid frequency for XTE J1550-564- Data points are la- 
beled with different marks according to the selection ex- 
plained in text (same symbols as in Fig. 7). 

ergy range and the PDS were fitted with Lorentzian mod- 
els plus a low frequency power-law. The PDS were rather 
complex: as centroid frequency we used the peak at the 
lowest frequency in the 1-20 Hz range. The 3-20 keV PCA 
energy spectra were fitted with a disk-blackbody plus a 
power law, although the presence of an excess around 6-8 
140 k e y was recognized. 





* - 
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QPO (Hz) 

Fig. 9. Plot of Power Law Photon Index versus QPO cen- 
troid frequency for XTE J 1748-288 (diamonds) and 4U 
1630-47 (filled circles). 

The points corresponding to this source are shown 
in Fig. ^J. It is possible to recognize a behaviour similar 
to those shown by GRS 1915+105 and XTE J1550-564, 
although in this case the correlation is not as smooth. 
The points corresponding to the end phase of the out- 
burst have a rather constant T around 1.7-1.8, while the 
QPO frequency decreases from 3 Hz down to 0.2 Hz (data 
from Tomsick & Kaaret 2000). Notice that the points cor- 
responding to the lowest frequencies (0.2 to 0.5 Hz) do 
match rather well the correlation followed by XTE J 1550- 



F. Vignarca ct al.: Tracing the power-law with the QPO frequency 



9 



564. At higher QPO frequencies there is a correlation be- 
tween the two quantities, with a rather large scatter and 
a possible indication of a turn-off. 

4. Summary and discussion 

Our results can be summarized as follows. For the \ class 
of GRS 1915+105, when the energy spectrum is rather 
hard and consists exclusively of the power-law component 
and the power density spectrum shows a clear 1-10 Hz 
QPO, there is a strong positive correlation between QPO 
ccntroid frequency and photon index of the power law. 
At the highest values of the QPO ccntroid frequencies, 
this correlation flattens and appears to start reversing. 
Published data from other sources indicate that this gen- 
eral correlation exists for a number of sources, and indeed 
in some cases there is a turnoff in the relation for high 
QPO frequencies. The correlation is not exactly the same 
for all sources and in the case of XTE J1550-564 there is 
evidence for different turnoffs during different intervals of 
the source outburst. 

The turnoff in the correlation is difficult to understand. 
If the QPO frequency is associated to a specific radius in 
the accretion disk, only the rising part of the correlation 
can be explained with a simple Comptonization model. 
If the QPO-radius relation also shows a reversal (see for 
instance Rodriguez et al. 2002 and Varniere et al. 2002), 
then a turnoff like that presented here can be produced. 

The source for which we have more data, and the 
largest range in QPO frequencies is XTE J1550-564. From 
Fig. 8 one can see that not only the points from the first 
two sections of the outburst (see Fig. 7) follow a very 
narrow correlation, but also that all subsequent points lie 
below this correlation. This suggests that the correlation 
followed in the first part, when the overall PDS was in the 
form of a flat-top noise component plus a QPO (Cui et 
al. 1999; Rcmillard et al. 2002b), is perturbed in a very 
specific way when the source changes state and a powerful 
disk component appears in the energy spectrum (Sobczak 
et al. 1999b, 2000a, b). There is increasing evidence that 
the power-law component observed in the energy spectra 
of BHCs in the Low/Hard state is physically different from 
that seen in the Very High State (see e.g. Zdziarski et al. 
2001). It is therefore possible that the turnoff is due to a 
gradual transition between these two regimes. 

It is interesting to discuss the position of the turnoff in 
the QPO-r relation. Figure 10, where all the sources are 
shown together, suggests that: (a) XTE J1550-564 and 
GRO J1655-40 have a similar turnoff frequency around 10 
Hz; (b) XTE J1748-288 has a higher turnoff frequency, if 
any, larger than 30 Hz; (c) GRS 1915+105 (see Fig. 5) has 
a lower turnoff frequency at about 5 Hz. It is tempting to 
compare these values with the relative values of the esti- 
mated dynamical masses for the compact objects in the 
systems. Excluding XTE J1748-288, for which no dynami- 
cal mass estimate is available, we see that GRS 1915+105 
has a mass roughly twice and a turnoff frequency half 
that of GRO J1655-40, as one could expect with a sim- 



ple mass scaling. However, the recent mass estimate for 
XTE J1550-564 is around 10© (Orosz ct al. 2002), from 
which we would expect a turnoff frequency lower by a fac- 
tor of ^2 with respect to that of GRO J1655-40, which is 
evidently not observed. 

In conclusion, we have found that in GRS 1915+105 
and a number of other BHCs, the centroid frequency of the 
low-frequency QPO and the photon index of the power- 
law component in the energy spectra follow a very similar 
kind of correlation. In particular, for some sources there is 
evidence for the presence of a turnoff in this correlation. 
If the QPO frequency is associated to a specific radius of 
the accretion disk, this provides an important clue to the 
possible origin of the power-law component. More data are 
needed in order to understand the physical nature of the 
turnoff, but these results stress the importance of such 
simple spectral/timing correlations for black hole candi- 
dates. 
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Fig. 10. Plot of Power Law Photon Index versus QPO centroid frequency. Plateau results from GRS 1915+105 are 
plotted in full circles; empty crosses are data from GRO J1655-40; empty circles are data from 4U 1630-47; diamonds 
are data from XTE J1748-288; empty triangles are data from XTE J1550-564. 



